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Abstract. We present new results on CO depletion in a sample of nearby pre-stellar cores, based on observations 
of the millimeter C 17 and C ls O lines and the 1.3 mm dust emission with the IRAM 30 m telescope. In most 
cases, the distribution of CO is much flatter than that of the dust, whereas other tracers, like N2H + , still probe 
the latter. In the centre of these objects, we estimate CO to be underabundant by a factor 4-15 depending on the 
cores. The CO underabundance is more pronounced in the central regions and appears to decrease with increasing 
distance from the core centre. This underabundance is most likely due to the freezing out of CO onto the dust 
grains in the cold, dense parts of the cores. We find evidence for an increase of the CO depletion degree with the 
core density. 
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1. Introduction 

Of all the molecules present in the interstellar medium, 
CO (and its isotopes) has been the object of particular 
attention. It has been widely used to study the mass dis- 
tribution and kinematics of star forming molecular clouds. 
Moreover, CO is thought to play a major role in the ther- 
mal c ooling of molecular clouds (e.g., Goldsmith & Langer, 
1978 ). In very dense and cold regions of the interstellar 
medium the CO molecules tend to leave the gas phase 
and condense out onto the dust grains, di rectly affecting 
the thermal balance of clouds (Goldsmith, 2001 ). CO also 
participates in a complex chemical network and its de- 
pletion leads to changes in gas composition and chemical 
reactions. These, in turn, can influence the contraction of 
molecular clouds through ambipolar diffusion (e.g. Shu et 
al. 1997 ), via the ionisation fraction (Caselli et al. 1995 ) or 
the g a s phas e deuteration of molecules (Roberts & Millar 
2000| a, |200C| b). 

This phenomenon affects in particular pre-stellar con- 
densations, which are characterized by very low tem- 



peratures (~ 10 K) and high central densities (~ 10 4 — 
10 5 cm -3 ). These objects, still starless, are thought to 
evolve slowly via ambipolar diffusion towards higher de- 
grees of condensation until they become unstable and col- 
lapse under their own gravity to form a protostar. This 
pre-collapse phase is believed to strongly influence the fol- 
lowing stages of the star f orma tion process - mass infall 
rate (e.g. Henriksen et al. 1997 ), fragm entation, mass of 
future star (e.g. Bacmann et al. 2000 ). It is therefore of 



prime importance to characterize the physical and chem- 
ical conditions in the gas and dust of pre-stellar conden- 
sations which have not yet started to collapse. The study 
of CO depletion is in this context particularly relevant 
due to the implications of e.g. modified thermal balance 
or contraction by ambipolar diffusion induced by a de- 
creased abundance of the molecule. CO depletion in dens e 
cores has alread y been discussed by Kramer et al. ( 1999 ) 
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Willac y et al. ( |1998j ), and Caselli et al. fll999j p002a , 
|2002b ) who focussed on L1544, a core similar to the ones 
we discuss in this work. Here, we studied the distribution 
and abundance of CO via its rarer isotopes C 17 O(l-0) and 
C 18 O(l-0), and of N 2 H+(l-0) in a sample of pre-stellar 
cores, using the 1.3 mm continuum as a tracer of the H 2 
column density. Our aim was to quantify the depletion 
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and to link it with physical properties of the cores. We 
first briefly present the observations in Section |[ then the 
spatial distributions and column density estimations in 
Section^, discuss our results and conclude in Section^. 

2. Observations and data reduction 

We have selected a sample of 7 pre-stellar cores previously 
studied by Bacmann et al. ( 200C| , cf. Table T he core 
L1544, studied in detail by Caselli et al. ( 1999| , |2002a , 



2002b), has also been included for comparison. 



Observations in the millimeter lines of N 2 H + (l-0), 
C 17 O(l-0), C 18 0(2-l) and (1-0) were made in August 
2000 and 2001 at the IRAM 30 m telescope. The spectra 
were taken in position switching mode with a 20" sam- 
pling along cuts in chosen directions passing through the 
maximum of the 1.3 mm emission. The reference positions 
were checked to be free of C 18 O(l-0) emission. We used 
an autocorrelator with a resolution of 20 kHz as back- 
end. The spectral resolution was degraded so as to pro- 
vide a kinematic resolution of ~ O.lkms -1 at all wave- 
lengths. Pointing was checked typically every 1.5 hour and 
found to be better than ~ 3". The data were reduced us- 
ing the CLASS package. The line fluxes are expressed in 
units of main-beam brightness temperatures. The uncer- 
tainties in the fluxes are f» 10%. The beam size (~ 22" at 
110 GHz) and th e efficiencies (FWHM) are taken from the 
IRAM web page http: / /www. iram.es/ . Previous mapping 



of the cores in the 1.3 millimeter continuum is described 
in Bacmann et al. ( 2000] ). 



3. Results 

3.1. Distribution of molecular and dust tracers 

We compared the spatial distribution of CO, N 2 H + and of 
the dust traced by the 1.3 mm continuum emission. To this 
effect, we plotted for each core the integrated intensity of 
the molecular lines (C 18 O(l-0), C 17 O(l-0) and N 2 H+(1- 
0)) together with the millimeter continuum flux, along di- 
rections passing through (or close to) the maximum of 
the 1.3 mm continuum map. To obtain the integrated line 
intensities, we first fitted a gaussian (or a series of gaus- 
sians in the case of hyperfine structures) to the considered 
line, and then evaluated the area under the gaussian. The 
results obtained are similar if we consider the integrated 
signal between two given velocities. The resolution of the 
1.3 mm continuum map was degraded to the resolution of 
the molecular line observations, 22". The resulting plots 
are presented in Fig. [I] for all cores except L1544. Whereas 
the continuum is sharply peaked (except for L1709A), the 
distribution of the C 18 O(l-0) and C 17 O(l-0) integrated 
intensities is flatter (e.g. L1689B and L310) or completely 
fails to trace the strong continuum peak (e.g. L429 and 
OphD). The core L328 is particular in the sense that the 
CO isotopes are also sharply peaked (their distribution is 
only slightly broader than that of the dust: ~ 15% broader 
at half- maximum) . For all cores, the distribution of N 2 H + , 



on the other hand, closely follows that of the dust traced 
by the continuum. 

These results tend to show that CO is underabundant 
with respects to H 2 (traced by the continuum) in the cen- 
tres of the cores. As Caselli et al. ( |l999| ) pointed out, this 
is most likely due to the depletion of CO onto the dust 
grains, because of the high densities and low temperatures 
in the inner parts of the cores. However, it can already be 
seen from the various behaviours of the plots such as those 
in Fig.|l] that the amplitude of this phenomenon, the de- 
pletion degree, is not the same in all the cores considered. 

3.2. Physical Conditions 

The average H 2 density at the emission peak in the cores 
was estimated to be of a few 10 5 cm~ 3 in the 11" telescope 
beam by Bacmann et al. (2000, see also Table |l|). Hence, 
the populations of the energy levels of CO and its isotopes 
are very close to thermalization. The ratio of the C 17 
(2-1) to (1-0) line integrated areas, calculated under the 
assumptions of optically thin lines and at LTE, therefore 
provides a good estimate of the kinetic temperature. We 
applied the same procedure to the C 18 lines under the 
same assumptions. We obtained consistent results between 
both isotopes, with values in the range 6-8 K for the C 17 
lines and in the range 7.5-10 K for C 18 0. The results are 
summarized in Table[l[ Such temperatures are similar to 
those found towards other pre-stellar and/or protostellar 
cores (Caselli et al. 1999; Fuller & Myers, 1993; Ladd, 
Fuller & Deane, 1998). 

We evaluated the CO column density using the rare 
isotope C 17 in order to make sure our results were not 
affected by opacity effects. The opacity of the C 18 O(l-0) 
transition was determ ined following the method described 
by Ladd et al. ( |1998|) . Suppos ing a [ 17 Q1/[ 18 Q1 ab undance 
ratio of 3.65 (Wilson & Rood |1994| ; Penzias |l98l| ) , we de- 
rived opacities lower than 0.1 for Oph D, L328 and L310, of 
1.1 for L429, 1.3 and 1.5 for L1709A and L1689B, respec- 
tively. According to the aforementioned C 18 O(l-0) opac- 
ities, the C 17 O(l-0) transition is optically thin for all of 
the cores. This validates our determination of the kinetic 
temperature from the C 17 lines. Although the opacity 
effects are small in the general case, we prefered to use 
the same tracer C 17 for all cores for more accurate com- 
parisons. The C 17 column density was determined under 
the same hypothesis and adopting as kinetic temperature 
the values derived previously. The values of Nqitq deter- 
mined at the continuum emission peak are presented in 
Table § 

The H 2 column density was determined from the dust 
thermal emission at 1.3 mm, which is largely optically 
thin. We degraded the resolution of all our continuum 
maps to the resolution at the frequency of the C 17 0(1- 
0) line, so as to compare H 2 and C 17 column densities 
in the same beam. The dust temperature was taken to 
be 7.5 K in all cores at their centre, as suggested by the 
recent works of Zucconi et al. (2001) and Evans et al. 
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Table 1. Comparison of C 17 and (calculated from dust) H2 properties. The H2 densities are taken from Table 3 of 
Bacmann et al. (2000). Errors on the temperature are typically ±1K. Uncertainties on Nn 2 and n(H2) are a factor of 
~ 2, those on X and / a factor of ~ 2.5. 
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(2001). This is slightly lower than in the envelope, where 
values of 11 K were determined from ISOPHOT obser- 
vations (Ward-Thompson et al., 2001). Following Andre 
et al. (1996), we have adopted a dust opacity coefficient 
Ki.3 mm = 0.005 cm 2 g _1 , appropriate for gas in the central 
pre-stellar regions, where the H2 density is ~ 10 5 cm _3 , 
though the act ual va lue is uncertain by a factor of 2 
(Preibisch et al. , |l993[ ) . The values of the H2 column den- 
sity at the maximum of the 1.3 mm continuum emission 
are shown in Table! 



4. Discussion 

The results presented above show strong evidence for CO 
depletion onto the dust grains in the densest parts of the 
molecular cloud cores. We present in Table |] the ratio X 
of the C 17 and the H2 column densities. For all cores, the 
abundance X found is lower than the "canonical" abun- 



dance determined by Frerking et al. ( 1982 ) towards dark 
cores X can — 4.8 x 10~ 8 and varies from 4.5 times smaller 
to nearly 16 times smaller than X can . The CO molecules 
thus appear undcrabundant with respect to H2 at the cen- 
tre of the pre-stellar cores. The values found suggest differ- 
ent depletion degrees for the different cores, as was already 
hinted at by the cuts of Fig.[fJ. The depletion factor /, de- 
fined as X can /X, is also shown in Table 1. For L1544, 
we derive a depletion factor of 14, the difference with 



the value of 10 found by Caselli et al. (1999) being due 
to our assumption of a colder central dust temperature. 
Adopting a higher value for the dust absorption coefficient 
Ki.3mm, representative of denser environments would de- 
crease the depletion factor to 6 in OphD, L429, L310. This 
case is extremely unlikely as it requires typical densities 
of - 10 6 cm" 3 over a size of 20" (3200 AU). The depiction 
factor would also be a factor of ~ 2 smaller if the central 
temperature is taken to be 11 K instead of 7.5 K. These 
results considerably extend the previous works on deple- 
tion in pre-stellar cores by Kramer et al. (1999), Caselli 
et al. (1999) and more recently by Tafalla et al. (2002). 
The latter report similar depletion factors for L1544 and 
4 other pre-stellar cores in the Taurus molecular cloud. 

The effect of CO underabundance is still detected at 
large scale, with a smaller amplitude, in regions covering 
large parts of the lower density gas, where the "standard" 



value of Ki.3 mm (0.005 cm 2 g _1 ) is a reasonable approxima- 
tion. For each core, we derived the C 17 and H2 column 
density along our cuts at 70% and 50% of the maximum of 
the continuum emission. The values of 70% and 50% were 
chosen as a compromise between having regions of lower 
density and not being affected by the loss of sensitivity of 
the millimeter co ntinuu m at the edges of the 1.3 mm maps 
(e.g. Andre et al. 1996 ). We assumed higher dust temper- 
atures of 10 K and 11 K to derive the depletion factors / 70 
and / 50 , respectively, since the outer layers are tho ught t o 
be warmer than the core centre (e.g. Zucconi et al. 2001 ). 
The abundance X 70 is 2 to 5 times smaller than X can , the 
CO underabundance is more pronounced in the densest, 
central regions. For all cores, the average depletion factor 
/ 70 is lower than /, and / 50 is lower than / 70 . 

Such variations can be easily understood if depletion 
is the cause of this underabundance of CO. Indeed, the 
less dense and warmer parts of the cores should be less 
affected by the freezing out of CO molecules onto dust 
grains. Within the objects, the variations in density are 
about one order of magnitude, based on the density fits 
of Tafalla et al. (2002), whereas the temperature varies by 
~ 50% at most. Hence, the variations of the depletion de- 
gree are probably dominated by the variations in density 
and not in temperature. We stress the fact that the rela- 
tive variations between /, / 70 and / 50 remain true even if 
the temperature is assumed uniform in the cores. This re- 
lation shows that the depletion factor increases with den- 
sity inside of the core, independently of any modelling of 
the internal core structure. The increase in / is particu- 
larly pronounced in the case of the densest cores L429 and 
OphD. 

We tried to relate the CO underabundance with the 
physical properties of the core. To this effect, we plot- 
ted the depletion factor / as a function of the mean den- 
sity in the core central region (Fig.§. The above plot 
indicates an increase of the depletion factor with den- 
sity, confirming the trend observed for individual cores. 
A linear regression of the form / cx n(H2) Q performed 
on the whole data set yields a ~ 0.1 — 0.5. The rather 
large range of values allowed for the exponent is mainly 
due to the point L310 which stands out rather markedly 
of the trend mentionned above. Indeed, leaving aside the 
L310 data, the points match rather well a relation of the 
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Fig. 1. Spatial distribution of the integrated line intensity of C 18 O(l-0) (filled triangles), C 17 O(l-0) (empty triangles), 
N 2 H + (l-0) (empty diamonds), and of the dust continuum emission (solid line) at 1.3mm, along cuts in R.A., Dec, or 
the diagonal Aa = AS (for L310 and L1709A) passing through the maximum of the continuum emission. 



form / oc rj(H2)°' 7 ~ ' 8 . We note that the observational 
uncertainties are much larger on L310 than on the rest 
of the dataset. First, the weather was noticeably unsta- 
ble when L310 was observed. Because the core is highly 
elliptical, the pointing fluctations induced by anomalous 
refraction could lead to underestimating the flux of the 
C 17 (1-0) line. Second, in deriving the average density 



tricity of the core. They took as size scale the dimension 
along its major axis, which leads to an underestimate in 
their derivation of the density. The L310 point should be 
actually shifted in the direction of a lower depletion de- 
gree and a higher density, as indicated by the arrow in the 
graph, so that the L310 point probably lies closer to the 
rest of the data. Therefore, the depletion degree appears 



in L310, Bacmann et al. (2000) neglected the high eccen 
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Density [10 5 cm" 3 ] 

Fig. 2. Variations of the depletion factor / = X can / X as a 
function of the mean central density. The parallellograms 
represent the error boxes for each point, accounting for 
the fact that the error on Ki.3 mm affects both density and 
depletion estimates simultaneously. The arrow indicates 
the direction the point for L310 would move with better 
density estimates and more reliable C 17 measurements. 
The dashed and dash-dotted lines represent the simple 
model / = exp(t/rd), where t& is the depletion timescale, 
for various core ages. 



to increase with density, following a power law close to 
/ oc n(H 2 ) a4 -°' 8 . 

However, one has to be cautious about the "exact" 
power law coefficient, given the uncertainties in the H2 
density, due for the most part to the uncertainty in Ki.3mm- 
The above relation was derived assuming a constant dust 
opacity coefficient and one may wonder how the exponent 
would be affected in the case fti.3 mm increases with den- 
sity. In such a case, the H2 column density differs from the 
evaluation by Bacmann et al. (2000), which, in turns, af- 
fects the estimation of the depletion degree too. We adopt 
a power law Ki.3 mm oc n(H2)' 3 , where ~ 0.2 under the 
standard assumption that Ki.3 mm increases by a factor of 
2 (from 0.005 to 0.01 cm 2 g _1 ) when the density increases 
from 10 5 to 10 6 cm" 3 (Preibisch et al. 1993). A simple al- 
gebra shows that the corrected slope of the / — n(H2) rela- 
tion is then a(l + 0) — instead of a (a is estimated, like 
before, supposing Ki.s mm constant). The corrected slope is 
less steep when a < 1, and it is positive if a is larger than 
0/(1+0) ~ 0.16. The slope a derived previously (0.4-0.8) 
lies well above this value of 0.16, therefore the corrected 
slope is positive (i.e. the depletion factor is found to in- 
crease with density, even in the case Ki.3 mm also increases 
with density). We conclude that the observed variations 
of the depletion factor with density cannot be accounted 
for by variations in Ki.3 mm only. Estimates of H2 column 
density usin g the Near Infrared Color Excess method (cf. 
Bcrgin et al. [2001 ) that are independent of Ki.3 mm should 
help confirm this trend. 

Assuming that a CO molecule in the gas phase has a 
probability equal to one to stick onto the grain after a col- 



lision (e.g. Leger 1983| ), the CO depletion is expected to 
vary roughly as exp (t/8000yr) exp(n(H 2 )), where rt(H 2 ) is 
expressed in units of 10 5 cm -3 . We plotted in Fig. || this re- 
lation for t = 2 x 10 3 yr and f = 2x 10 4 yr (dashed lines), 
as well as a "fit" for t = 8 x 10 3 yr (dash-dotted line). 
The observed scatter between the data and the model 
may well reflect the difference in evolutionary stages of 
the pre-stellar cores studied here. The dash-dotted model 



is consistent with the results of Caselli et al. (1999) on 
L1544. We note that according to the models, the core ages 
are found relatively small compared to ambipolar diffusion 



timescales (Ciolek & Basu |2001| , Caselli et al. |2002b| ). We 
stress however that this simple model docs not take des- 
orption from the grains into account, and therefore the 
ages might be underestimated. 

A more detailed study of the relation with depletion 
and density within the cores will be published elsewhere. 
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